Polymer chain generation for coarse-grained models using radical-like polymerization 
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An innovative method is proposed to generate configurations of coarse grained models for polymer 
melts. This method, largely inspired by chemical "radical polymerization", is divided in three stages: 
(i) nucleation of radicals (reacting molecules caching monomers); (ii) growth of chains within a 
solvent of monomers; (iii) termination: annihilation of radicals and removal of residual monomers. 
The main interest of this method is that relaxation is performed as chains are generated. Pure mono 
and poly-disperse polymers melts are generated and compared to the configurations generated by 
the Push Off method from Auhl et al. [J. A detailed study of the static properties (gyration radius, 
mean square internal distance, entanglement length) confirms that the radical-like polymerization 
technics is suitable to generate equilibrated melts. The method is flexible, and can be adapted to 
generate nano-structured polymers, namely diblock and triblock copolymers. 

PACS numbers: 61.41. +e Polymers, elastomers, plastics; 82.20.Wt Computational modeling, simulations; 
82.35 Jk Copolymers, phase transitions, structure; 82.35.Lr Physical properties of polymers. 



I. INTRODUCTION 

Molecular simulation is becoming an increasingly pop- 
ular tool for the investigation of mechanical and thermo- 
mechanical properties of polymer materials. It can be ap- 
plied to investigate the properties of homopolymer sys- 
tems as well as to nanostructurcd copolymers or polymer 
based nanocomposites, and to gain a microscopic under- 
standing of the properties of these technologically impor- 
tant materials. 

The main issue is to understand relations between 
polymer nanostructure and, in particular, mechanical 
properties. In order to bridge the gap between micro and 
macro scales, coarse grained molecular dynamics, where 
each "bead" represents several monomers, are becoming 
a standard tool. They allow for an investigation of qual- 
itative and quantitative issues not directly accessible to 
experiments, while remaining affordable in terms of com- 
putational costs. 

Investigating structure-property relations in polymeric 
systems, however, requires the preparation of equili- 
brated melts with long and entangled chains. Above the 
glass transition, equilibrium can, in principle, be achieved 
using long Molecular Dynamics (MD) or Monte Carlo 
(MC) simulations. The situation gets difficult for long 
chains with relaxation times that can soon exceed the 
typical simulation duration of a few nanoseconds, and for 
nanostructured polymers (e.g. block copolymers), where 
equilibration times, even for short chains, are too long to 
use MD or MC to produce equilibrated melts. 

For long polymer chains, hybrid methods combining 
MD and MC inparticular the so called "double bridg- 
ing" algorithm [2j , have been used to generate well equi- 
librated melts. These algorithms, apart from their tech- 
nical complexity, are not particularly well suited for ex- 
tension to more complex architectures. 

The objective of our contribution is to propose an inno- 
vative method for polymer chain generation, (i) based on 
a realistic approach close to radical polymerisation^, 
(ii) particularly adapted to generate non linear architec- 



tures (branched polymers, star polymers, copolymers,...) 
and/or polydispcrsc chains; (iii) providing equilibrated 
melts. 

This method, called "radical-like polymerization" will 
be tested on different system types (mono- and poly- 
disperse homopolymers). It will be also compared to more 
classical Push Off methods [!, Q , which are based on a 
two steps process: (i) random gaussian chain generation 
and (ii) equilibration. Systems resulting from step (i) are 
usually quite far from equilibrium as chains interactions 
are not taken into account, requiring thus long equilibra- 
tion times (step (ii)). 

The main idea of radical-like polymerisation, is that 
chains are partially relaxed simultaneously while poly- 
merization is achieved. 

The manuscript is organized as follows. Section |TT] de- 
scribes the method. In section [TTTT we apply the method 
to several types of homopolymer melts, and show how 
it can be tuned to obtained well equilibrated melts at a 
relatively low computational cost. Finally, we point out 
that the radical-like polymerization method is suitable for 
simulating block copolymers, and give a preliminary il- 
lustration of this in section [TVl 



II. 



DESCRIPTION OF SYSTEMS, 
METHODOLOGY 



AND 



Our simulations are carried out for a well established 
coarse-grained model [5| in which the polymer is treated 
as a chain of N = J2 a beads (where a denotes 
the species for block copolymers), which we refer to as 
monomers, of mass m = 1 connected by a spring to form a 
linear chain. The beads interact with a classical Lennard- 
Joncs interaction: 



TJ Q/3 (r) = I 4eaf} ( f7Q / 3 / r ) 12 ~ ( a af}/r) 6 ,r<r c 
I , r>r c 



(1) 
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where the cutoff distance r c = 2.5a a ^. In addition to |T]), 
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Parameters Signification 



0.5fci?gln(l- (r/Ro)' 



r<R 
r>R 



monom 



U F ene(t) 



The parameters are identical to those given in Ref. 
namely k — 30e Q ^/cr^ and Rq = l.§cr a p, chosen so 
that unphysical bond crossings and chain breaking are 
avoided. All quantities will be expressed in terms of 
length (Tii = 0") energy e aa = e and time tlj = ma 2 /e. 

Newton's equations of motion are integrated with 
velocity- Verlet method and a time step 5t = 0.006. 
Periodic simulation cells of cubic size L containing M 
chains of size N where used under a Nose-Hoover baro- 
stat, i.e. in the NPT ensemble. The pressure is fixed to 
P = 0.5e/a 3 



A. Radical-like polymerization 

1. Algorithm 

The radical-like polymerization method is directly in- 
spired by the classical polymerization phenomenon with 
a protocol based on three stages: 

• starting: a radical (active molecule that interacts 
with monomers) is created by an active molecule 
A (A — > P*) and interacts with a first monomer 
P* +M -> PM*, 

• propagation: the radical captures a new monomer 
and moves to the chain end PM* + M -> PMM* 

• termination: four main mechanisms of termina- 
tion can usually be identified in polymerization 
reactors: (i) two radicals can annihilate leading 
to two separated polymer chains (PM . . . M* + 
PM . . . M* -> PM ...M + PM . . . M) (dispro- 
portination); (ii) two radicals can annihilate leading 
to one polymer chain (PM ...M* + PM . . . M* — > 
PM ...MM ...MP) (coupling); (iii) a radical can 
be transferred to another monomer leading to a new 
growing chain (trans fert) or annihilated by some 
defect. Radicals can also remain active and chain 
growth is stopped only when all monomers have 
been consumed, as in (living polymerization). 

The radical-like polymerization process takes place in 
a solvent which is represented in our simulations as a 
Lennard-Jones liquid of N monom monomers. 

Note that the aim of our method is not to model in 
detail the polymerization process, but rather to take in- 
spiration from it. As a reminder, we give in Tab. [I] a sum- 
mary of relevant parameters fully describing the radical- 
like polymerization algorithm. 

The radical-like algorithm is then divided in four 
stages: 

1. Nucleation: each monomer has a probability p 
to be randomly functionalized as a radical. This 
probability p controls the number of chains M = 
p x N rnonom that will eventually be created. 



A, 
M 
N, 
N 



ngrowth 
UMDSbg 
nrelax 

ntot 



Total number of beads in the simulation box 

Total number of chains 

Final length for a chain i 

Desired chain length for isodisperse systems 

Nucleation probability 

Number of growth steps 

Amount of MD steps between each growth step 
Number of MD steps during relaxation phase 
Total number of MD steps 



Tab I: Relevant parameters used in the radical-like polymer- 
ization algorithm. 



amount of growth steps n grow th, defined initially, 
controls the maximum chain length N^max = 
rtgrowth- Note that this allows us, as mentioned pre- 
viously, to mimic the polydispcrsity associated with 
living polymerisation. This stage of the process is 
schematically depicted on Fig. [TJ 

3. Relaxation: this is an essential ingredient of the 
method. Between two successive growth steps, a 
radical is allowed to explore its neighborhood dur- 
ing riMDSbg MD steps. This is equivalent to let a 
chain evolve in the solvent and explore a part of its 
conformational phase space in situ while polymer- 
ization is taking place, hence permitting a partial 
relaxation. 

4. Termination: For polydisperse systems the gen- 
eration procedure is stopped after a fixed number 
of growth steps n growt h- To produce a monodis- 
perse system, the process is stopped only when 
each chain has reached a desired size N , whatever 
the number of growth steps. Naturally, the time 
elapsed before termination will depend on the ratio 
(N x M)/N monom which we took near 80%. Finally 
the residual monomers (or solvent) are removed and 
the system is relaxed to reach at constant pressure. 




Fig 1: (Color online). Growth step during the radical-like 
polymerization algorithm. A radical (white) is randomly as- 
signed one of its first monomers neighbors (blue ones, num- 
bered from 1 to 4) to create a new covalent bond and increase 
the local chain length Ni. 



2. Growth: radical (index i) randomly chooses one 
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be included trough a preferential choice of the neighbors 
during the growth step. Three types of systems were gen- 
erated using the radical-like polymerization process: 

• Cold: pure polydisperse melt. The polymeriza- 
tion procedure involve a finite value n growt h of 
growth steps but without coupling the system to 
a heat bath by imposing riMDSbg = 0, thus pre- 
venting any relaxation between growth steps. 

• Hot: pure polydisperse melt. The number n growt h 
of growth steps is also fixed to a finite value, but for 
this kind of polymerization, the system is coupled 
to a heat bath by fixing a finite number riMDSbg of 
relaxation steps between each growth step, and set- 
ting MD parameters using a Nose-Hoover barostat 
with ksT = 2 and P = 0.5. For this kind of pro- 
cedure, the polymerization process is stopped once 
the number of growth steps is reached. 

• HotMono: pure isodisperse melt. For this kind 
of process, the number of growth steps is infi- 
nite. Practically, growth stage occur until all chains 
reach the desired size N. The system is coupled to 
a heat bath (Nose-Hoover thermostat and barostat 
with fcfiT = 2 and P = 0.5) during the relaxation 
stage. riMDSbg MD steps arc performed after each 
growth step. 

In the next section, these three types of generation 
processes will be tested and compared. The monodisperse 
HotMono generation procedure will also be compared to 
the more classical Push-Off techniques 0, [f|. 

Within the Push- Off framework, chains are gener- 
ated randomly in the simulation box without consider- 
ing excluded volume Q. Thus, Lennard- Jones interac- 
tions for non-bonded monomers cannot be introduced 
immediately because chains spatially overlap. To bypass 
this difficulty, modified Lennard- Jones potentials (Slow 
Push Off), or intermediate soft repulsive potential (Fast 
Push Off) arc then introduced, and eventually replaced 
by the LJ potential. Due to its relative simplicity, this 
method has been widely used in the literature to gener- 
ate monodisperse systems. 

We refer to [l[ for details and discussions about FPO 
technics. In our implementation, the systems generated 
with FPO (M = 200 chains with chain length of N = 
200) are relaxed during 10 7 MDS for systems o un- 
der Nose-Hoover thermostat (ksT = 2.0) and barostat 
(P = 0.5). It has to be noticed that this quite easy pro- 
cedure is known to create significant distortions in the 
chain statistics on length scales comparable to the tube 
diameter [H,@,0]; requiring thus relatively long relaxation 
times. Consequently chain length is generally limited to 
N < 400. 



2. parameters 

The values of the parameters used in our generation 
processes and subsequent simulation for the three types 
of protocols are summarized in Tab. [TTJ For polydisperse 
systems, e.g. Cold and Hot, the min and max values of 
the chain length distribution are also quoted in the same 
table, and will be discussed below. 



melt. This parameter can be considered as a control pa- 
rameter for the exploration of configurational phase space 
during growth, at a given temperature and pressure. 





^growth 


n-MDSbg 


M 


N 


(JV)(t->oo) 


Cold 


350 





184 


[50; 344] 


172 


Hot 


6.7 10 4 


10 


215 


[56; 390] 


226 


HotMono 


10 5 


10 


215 


200 


200 




10 5 


300 


215 


200 


200 



Tab II: Parameters used to simulate the different radical-like 
polymerization processes discussed in text, during the gener- 
ation stage. 



3. Structural characterization 

Three types structural parameters have been investi- 
gated to control the state of relaxation of polymer melts: 

• the mean gyration radius (r g ) defined by: 

M'-t ^'-™''" (3) 

i=l 

where r] is the position of the j-th atom of the i-th 
chain, (r 1 ) is the center of mass of chain i and N l 
is the size of chain i. 

• the Mean Square Internal Distance (MSID) (r 2 )(n) 
is the average squared distance between monomers 
j and j + n of the same chain. It is defined by: 

(r>>=s Egrw -lr )2/(w '' n> w 

1=1 

Note that the MSID (r 2 )(n) is a function of n and 
(b 2 ) 1/2 = v /(r 2 )(l) is the mean bond length. 

• the primitive path: the Primitive Path Analysis is 
a powerful tool to investigate the distance between 
chains entanglements. It is a key parameter, that 
controls the mechanical or rhcological properties of 
the polymer melt. Section fill CI will be devoted to 
the PPA. 

III. RESULTS FOR A HOMOPOLYMER MELT 

A. Dynamics of the polymerization process 

A preliminary study is devoted to the growth dynam- 
ics of polydisperse system, namely Hot and Cold meth- 
ods. In Fig. [2l the mean chain length (N) is plotted as a 
function of the number of growth steps preformed during 
polymerization. It is worth noting that polydispersity has 
spontaneously appeared as a result of the growth process. 
We observe that both methods display the same evolu- 
tion: a rapid increase followed by a saturation due to the 
lack of available monomers. However, the Cold procedure 
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Fig 2: (Color online) Mean chain length size (N) (symbols) 
and polydispersity index I p (lines) evolution during polymer- 
ization stage versus the number of growth step, and for the 
two Hot and Cold simulated isodisperse systems. Also plot- 
ted is the standard deviation crjv represented by vertical bars 
centered on symbols. Inset: Size distribution P(N) for the 
same systems at the end of the generation procedure. 



Fig 3: (Color online) Generation stage: evolution of the mean 
gyration radius (r g (t)) normalized by the average bond length 
(b 2 (t)) 1/2 and averaged over all chains. Generation exhibits 
two distinct stages: (i) a pure growth stage characterized by 
a growth kinetics; (ii) a saturation stage where gyration 
radii reach a plateau value. A value of UMosbg = 10 has been 
used for Hot and HotMono methods (see Tab. Hi) . 



The standard deviation = \J (N 2 ) — (N) 2 is also 
indicated for both systems with vertical bars centered on 
respective symbols. 

The final size distributions at the end of the generation 
procedure P(N) are plotted for both systems on inset of 
Fig. [5] As expected, the peak is shifted towards larger 
sizes and is slightly narrower for the Hot system. 

In our simulations the polydispersity index I p = 
Mw/Mn 1 , is accessible through the ratio I p = 
(N 2 ) I (N) 2 . The final polydispersity index is a little lower 
for the Hot system (around 1.057) than for the cold sys- 
tem (around 1.103). Again, this is probably due to ther- 
mal mobility which allows smaller chains to find new 
monomers to continue the growth. 

Our generation procedure, which is very close to liv- 
ing polymerisation (see section Hi Al) leads to polydisper- 
sity indexes that are reasonably close to experimental 
ones resulting from living polymerization (typically of or- 
der 1.3), which gives us confidence in the physical back- 
ground our the radical radical-like polymerization algo- 
rithm. Moreover, it would be very easy to slightly modify 
our method to simulate other kind of polymerization pro- 
cesses which would lead to higher polydispersity. Experi- 
mental values of polydispersity index can reach a value of 
10 or more for classical polymers where coupling, trans- 
fert or disproportination are involved (see section III Ap . 

In order to quantify the evolution of the structural 
properties of chains during production runs for the Hot, 
Cold and HotMono methods, we also investigated the 
evolution of the mean gyration radius (r g (t)) normalized 
by the mean bond distance (b 2 ^)) 1 / 2 during the growth 
(Fig. [3]) and relaxation (Fig. [3J stages. Such evolutions 
are investigated for the three systems (riMDSbg = 10 
for Hot and HotMono during the generation stage - see 
Tab. OH). 

In figure [31 we observe that the generation proceeds 
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Fig 4: (Color online) Relaxation stage (e.g. after polymeriza- 
tion): evolution of the mean gyration radius as a function of 
the number of MD steps necessary to reach a total number 
n to t = Ugrowth x nuDSbg + n re i ax = 10 7 MD steps. Fast Push 
Off (FPO) and HotMono methods converge to the same value. 



in two distinct stages: (i) a pure growth stage charac- 
terized by a growth kinetics; (ii) a saturation stage 
where gyration radii reach a plateau value. The power 
law simply means that during stage (i), each growth 
step is successful and lead to an increase of the chain 
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As r g cx N 2 , we obviously get 
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In figure 01 the time evolution of the mean radius 
of gyration for the Cold, Hot, HotMono and also Fast 
Push Off (FPO) are compared during the relaxation 
stage. The radius of gyration is plotted versus the num- 
ber of MD steps necessary to reach a total number 
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nal values of gyration radii depend on mean chain length 
N: the Cold method, which gives the smallest final mean 
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Fig 5: (Color online) Growth and relaxation stages: evolution 
of the mean gyration radius as a function of the mean chain 
size during growth (curves) and relaxation (vertical arrows) 
stages. Data from Kremer[5| are also represented. They pre- 
dict a Nft) 1 / 2 dependance. After the removal of remaining 
monomers and 10 7 MD relaxation steps, all generation tech- 
nics are in very good agreement with Kremer's results. 



radius. Finally, the Hot method, which gives the largest 
final mean chain length (N = 226) leads to the largest 
mean gyration radius. 

In order to investigate the evolution of the chain size 
as a function of chain length during the growth and re- 
laxation stages for all polymerization methods, we plot- 
ted (r g (t))/(b 2 (t))^ 2 versus (N) on Fig. [SJ In this fig- 
ure, relaxation process (at constant N) is represented by 
vertical arrows. We also plotted in this figure data from 
KremerJEJ resulting from long time equilibration, which 
predict a N 1 / 2 dependance. 

After the removal of remaining monomers and 10 7 MD 
relaxation steps, all generation methods [Cold, Hot and 
HotMono) are in very good agreement with Kremer's tar- 
get function (r g (t)) / (b 2 ^)} 1 / 2 versus (N). 

However, with Cold, Hot or HotMono) (with 
UMDSbG = 10), it seems that relatively long relaxation 
times (up to 10 7 MD steps) are necessary to reach Kre- 
mer's target function. Therefore, in what follows, the ef- 
fect of the number of MD steps between each growth step 
{nMDSbo) will be investigated. 

In Fig. [?51 the mean normalized radius of gyration is 
plotted versus simulation time for the generation of M = 
215 chains of expected length N = 200 at fc^T = 2 and 
P = 0.5. Two different values of nMDSbg are investigated: 
n-MDSbg = 10 and riMDSbg = 300. It can be observed 
that a larger value of UMDSbG slows down the growth 
kinetics, but leads to better equilibrated systems once 
growth is completed. For riMDSbG — 300, no equilibration 
MD steps are required to reach the radius of gyration 
obtained with the Fast Push Off (FPO) method. 

This shows that the chains generated here reach their 
equilibrium structure more rapidly for the protocole that 
spends more time during the growth stage 2 , thus pointing 
out the main interest of this algorithm: i.e. equilibration 
is occurring during generation, provided an appropriate 



Fig 6: (Color online) Evolution of the mean radius of gyra- 
tion as a function of time (in MD steps) during growth and 
relaxation stages: generation of M = 215 chains of expected 
length N = 200 at ksT = 2 and P — 0.5. Two different values 
of UMDSbG (the number of MD steps between each growth 
step) are compared. A larger value of UMDSbG slows down 
the growth kinetics, but leads to better equilibrated systems 
once growth is completed. For UMDSbG = 300, no equilibra- 
tion stage is required to reach the mean radius of gyration 
obtained with the Fast Push Off (FPO) method. 

compromise for the number of MD steps between growth 
steps is chosen. 



B. Comparison of chains structure for HotMono 
and FPO methods 

The structure of a polymer melt can be characterized 
by a wide variety of static or dynamic interchain and in- 
trachain correlation functions [H, H, H, 0, [HI EH which are 
more or less sensitive to the artifacts introduced by the 
preparation procedure and which equilibrate on different 
time scales. One may note that for fully flexible chains 
simulated in our model (only FENE + LJ interactions), 
i.e. without torsional barrier and bending stiffness poten- 
tials, the local monomer packing relaxes quickly, while 
deviations of chain conformations on large scale require 
large times to equilibrate. 

To validate our generation methods according to more 
"classical" techniques, we will be interested, in the follow- 
ing, by a measure of internal chain conformation, namely 
the Mean-Square Internal Distance (MSID) (r 2 )(n). This 
function, defined in Eq. (j4|) above, gives the average 
squared distance between two monomers belonging to the 
same chain, and separated by a subchain of n monomers. 

The MSID is shown in Fig. [7] for all simulated systems 
after the total number of MD steps n to t = n growt h x 
nMDSb g + n re i ax = 10 7 MD steps. Cold, Hot, HotMono 
(with nMDSbg — 10), and also Fast Push Off (FPO) seem 
to converge to the same configuration since they all fit 
nicely with the "target function" defined by Auhl[l[ as 
the signature of well equilibrated melts. 

Error bars in Fig. [7] arc estimated using the standard 
error function that includes the number of subset events 
taken into account to compute the MSID. As n reaches 
chain length A^ (n — > N) , less and less pairs of monomers 
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Fig 7: (Color online) Mean square internal distance (MSID) of 
generated melts measured after from long MD runs (10 7 MD 
steps). All the simulated systems (Cold, Hot, HotMono and 
FPO are compared to the target function of Auhl[l|. Error 
bars are calculated using standard error function on statistical 
samples. All methods lead to well equilibrated melts. 
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Fig 8: (Color online) Mean square internal distance (MSID) 
of mono-disperse melts. The effect of the number of MD steps 
between each growth step is studied. A larger value of riMDSbG 
leads to better equilibrated systems, whom which MSID fits 
nicely with FPO and the target function of Auhl[l(. 



for large n are not statistically significant. 

Once again, all our generation methods ( Cold, Hot and 
HotMono) lead to well equilibrated melts (according to 
the MSID criterion) after (i) generation, (ii) removal of 
remaining monomers and (iii) 10 7 MD relaxation steps. 

As far as the radius of gyration is concerned, we showed 
that the number of MD steps between successive growth 
steps (nMDSbc) has an effect on the final structure of the 
melt. Indeed, the number (riMDSbg) of relaxation steps 
between growth steps can be view as a relaxation process 
for chains during the polymerization stage. Therefore, the 
MSID of mono-disperse melts has been investigated for 
riMDSbg = 10 and n M DSb g = 300. 

On Fig. [5J MSID resulting from HotMono generation 
(with nMDSbg = 10 and nMDSbg = 300) are compared 
with MSID resultingfrom FPO generation and the tar- 



tively low equilibration time, it can be observed in Fig. [8J 
that the HotMono generation method with 300 MD steps 
between each growth step leads to relatively well equili- 
brated systems, even slightly better than FPO method. 
This corroborates previous results from Fig.[6l and points 
out, once again, the main interest of this radical-like gen- 
eration method: relaxation takes place while generation 
is performed. 



C. Primitive Path Analysis 

Entanglements between chains are an important topo- 
logical feature, that controls many dynamical properties 
of polymer melts. A practical tool for characterizing en- 
tanglements is the Primitive Path Analysis (PPA) , which 
will be the object of this section. 

Proposed by Everaers 12| with the aim of constructing 
a real space representation of dc Gcnnes' tube model, the 
PPA technique is an interesting tool for obtaining infor- 
mations about the density of entanglements which has 
not been accessible through other theoretical or direct 
experimental measurements. 

Recently Hoy and Robbins 0] applied this technique to 
quantify the effect of the generation procedure and of the 
relaxation for two types of generation methods, namely 
the FPO system and the Double-Bridging[l| relaxation 
technique. Following their idea, we apply this to our dif- 
ferent radical-like generation methods, first focusing on 
the comparison between HotMono and FPO method. 

The principle of PPA is the following: 



(i) We start with any given configuration, during the 
growth or in the final state, after or before the re- 
laxation step. 

(ii) The two chains ends are kept fixed, while the intra- 
chain pair interaction (covalent bonds) are shifted 
to get their minimum energy at a zero distance 
while increasing the bond tension in Eq. @ to 
k = 100; 

(iii) To prevent chain crossing [l3j|, weak bonds lengths 
have been monitored and limited to 1.2a. 

(iv) The system is then equilibrated using a Conjugate 
Gradient algorithm in order to minimize its poten- 
tial energy and reach a local minimum. 

(v) The contour length of the primitive path L pp is then 
the total length of the chain (the sum of all straight 
primitive path segments length). 

If no entanglement exists between chains, L pp should 
be equal to their end-to-end distance r ee . The presence 
of entanglements leads to L pp > r ee with a typical Kuhn 
length a pp = {r1 e )/L pp and an average bond length 
(b PP ) = L pp /N. The number of monomers in straight 
primitive path segments is then given by: 



Npp(N) 



Tl 

VP 



(5) 



For short chains without any entanglements, the prim- 
itive path length equal end-to-end distance leading to 
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Fig 9: (Color online) Evolution of the number N e (t) of 
monomers in straight primitive path segments along simu- 
lation times for isodisperse systems HotMono and FPO. Cal- 
culation were performed both during the generation (poly- 
disperse) phase and the relaxation (isodisperse) phase sep- 
arated by the vertical dashed line. Units of time are in r 
units, i.e. n to t X St. The horizontal line gives value for iV e 
from Sukumaran[l3j . 



several entanglements per chains exist, and N pp (N) = 
N e . 

The PPA analysis has been performed at different sim- 
ulation times (during generation and relaxation stages) 
and results are shown in figures [5] and 1101 

Fig. [5] displays the number of monomers in straight 
primitive path segments N pp = N e for FPO and Hot- 
Mono (riMDSbG = 10 and riMDSbg = 300) generation 
methods. The vertical dashed line separates the genera- 
tion and growth regimes. The horizontal line is the en- 
tanglement length N e from Sukumaran[l3j|, that is in ex- 
cellent agreement with our data. This asymptotic value 
is even reached during the generation stage for the Hot- 
Mono technique with nMDSbg = 300: the relaxation stage 
is not required for this system! 

The radical-like method appears to perform perform 
particularly well at the relaxation stage: the entangle- 
ment length do not deviate much from the asymptotic 
value for the HotMono system in comparison to the FPO 
method. 

The PPA analysis has also been implemented for poly- 
disperse Hot and Cold systems. Fig. [TU] shows the ratio 
N e (t)/(N(t)) for polydisperse systems (Cold and Hot) 
against simulation time. During the generation stage, the 
time scale is given in n growt h steps units, whereas given 
in n re i ax MD steps during the relaxation stage. 

For Cold system, generation/relaxation transition is 
represented by a dashed vertical line, while a dot-dashed 
line is used for Hot system. 

The same indicative value for the entanglement length 
N e /N from Sukumaran[l3| for chain length of size N = 
200 is also shown, and must be considered as a mean 
value for both polydisperse systems. 

Indeed, the mean chain length at the end of the gen- 
eration phase for Cold system is (N)cold(t — > oo) = 172, 
while for the Hot (N) Hot (t -> oo) = 226 (see table HE]). 

For the Cold method, the investigated ratio is almost 



Fig 10: (Color online) Ratio N e (t)/(N(t)) for polydisperse 
systems Cold and Hot against simulation time. Dashed (mid- 
dashed) vertical line separates generation to relaxation stages 
for the Cold (Hot) method. Also shown is the same ratio from 
Sukumaran[13| for chains length N = 200 as an indicative 
value. 

plex behavior. First, a power law decrease, as noted by 
dotted (oc i -1 / 2 ) and dotted-dashed line (oc t^ 1 ), is ob- 
served, until n growt h ~ 700, corresponding to a ratio 
N e (t)/(N(t)) - 1/3 nearly equal to results from [l]| for 
N = 200 homopolymer chain melts. 

In this regime, (N(t)) grows more rapidly than N e (t), 
and the growth process of each chain interacts with a 
stochastic background associated with the ensemble of 
growing chains. Thus, in this Rouse-like regime, topolog- 
ical constraints do not play a significant role and one may 
expect that chains with average length (N(t)) < N e ~ 
N/3 dominate the polymerization, following a Rouse- like 
chain dynamics. 

Following this regime, while (N(t)) still grows, a sta- 
bilization of the same ratio is observed. In this regime, 
N e (t)/(N(t)) < 1/3, and a slowing down is observed 
during chain growth dynamics. This new reptation-likc 
regime, corresponds to a dynamics where the surrounding 
medium topology limits transverse chains displacements 
around their own contour length. Chains with mean size 
(N(t)) > N e ~ N/3 follow this reptation-like dynamics, 
and the polymerization process is slowed down. While 
the longest chains are still growing, the average entan- 
glement length does not vary drastically, as one can see 
once the generation stage is finished, where the ratio 
N e (t)/(N{t))^N e /N. 

From all these results, it appears that our approach is 
validated as a method for generating equilibrated config- 
urations of homopolymer melts. In the following section, 
the radical-like algorithm will be used to generate block 
copolymers in a lamellar configuration. 



IV. APPLICATION TO COPOLYMER 
GENERATION 

In this section, generation of block copolymers will 
be investigated and our radical-like polymerization al- 
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from a Lennard- Jones liquid of monomers: 

1. Nucleation: Each monomer i has a probability p 
to be a radical of type A if, say, > L z /2 and B 
otherwise. 

2. Growth: As long as the chain does not reach the 
size N/2 (N(t) < N/2), growth is performed as in 
a homopolymer with a supplementary condition: 
addition of a new monomer j is possible only if it 
lies in the same region (zj > L z /2 for A chains and 
Zj < L z /2 for B chains). The interface situated at 
z = L z /2 is then impermeable: no chain can cross 
it. 

3. From one region to the other: Once a chain 
reach the critical size N/2 (N(t) = N/2, the growth 
within a lamella is stopped. A force is applied to at- 
tract the chain ends to the interface, and the condi- 
tion above is reversed: addition of a new monomer 
j is possible only if it lies in the opposite region 
(zj < L z /2 for A chains and Zj > L z /2 for B 
chains). Under this new condition, and once a radi- 
cal combines with a new monomer in the opposite 
region, it turns into the opposite species (A rad- 
ical becomes B radical and B radical becomes A 
radical). For chains, with N(t) > N/2, the growth 
is then continued with the impermeable interface 
condition: addition of a new monomer j is possible 
only if it lies in the same region (zj > L z /2 for A 
chains and zj < L z /2 for B chains). Growth a a 
chain occurs until its length reach the size N. 

4. Relaxation: As for homopolymers, a number 
n-MDSbg of MD steps is performed between each 
growth step, during which the systems is coupled 
to the heat bath at ksT = 2e aa and P = 0.5. 



Fig 11: (Color online) Snapshots of symmetric di-block 
copolymers generated using the radical-like co-polymerization 
method for M — 215 chains of length N = 200, and for two 
values of the Flory-Huggins compatibility parameter \N = 
396 (segregation regime - up), and \N — 4 (one phase re- 
gion: mixing in progress - down). 



viously assembled copolymers, Srinivas F| managed to 
obtain large scale demixtion in biological systems (self- 
assembled copolymers in water). Zang|l5[ used full- 
atomistic simulations based on dynamics density func- 
tional theory but their approach is limited to small sys- 
tem sizes. May be more adapted to block copolymer 
generation, semi-particle based methods such as Single 
Chain in Mean-held 0, [l?], [H, EH seem to be promis- 
ing. 

In the following, we present an alternative method 
based on an adaptation of the radical-like method to the 
particular case of symetric AB diblock where Na = Nb 
and N = Na + Nb- L x , L y and L z are the box sizes along 



anrl v Hi rnr-i i nnc 



With this procedure, all chains have the same length 
N = Na + Nb and Na = Nb- Systems are then relaxed 
at k B T = 0.5e aa and P = 0.5 during 10 6 MD steps. 

Values for excluded volume potentials Q] and [3] have 
been chosen as, €aa = £bb = 1-0 and the interaction 
with the solvent fixed to e as — 1.0 where a € [A, B] and 
s being the solvent molecular type. All a a p = 1.0 while 
potentials are truncated and shifted at r c = 2.ba a p. 

The order-disorder transition temperature is gov- 
erned by the product of %7V, where x = ( e AA + 
e bb — 2cab) I i^ksT) is the Flory-Huggins temperature- 
dependant interaction parameter characterizing the AB 
incompatibility. Symmetric diblock copolymers are ho- 
mogeneous at small xN value, but strongly heteroge- 
neous with ordered structure when xN exceeds, in mean- 
field theory, the critical order-disorder transition value 
xNodt ~ 10. Hence, as a first application of our radical- 
like copolymerization algorithm, we simulated such di- 
block copolymers in the two limiting case of weak segre- 
gation limit with xN = 4 ((eab — 0.99), and the strong 
one, with X N = 396 {e AB = 0.01). 

To observe box dilatation and inter-lamcllac distance 
relaxation, an anisotropic Nose-Hoover barostat has been 
used during 5.10 6 MD steps, in such a way that P x = 
P y = P z = 0, while the temperature was fixed to ksT = 
0.5e ace . 

o _r i_ a j.: 
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the upper frame of Fig. [TT], one can observe that the in- 
terface separating the two blocks is well defined and sta- 
ble, as expected in the strong segregation limit. On the 
contrary, in the right panel, the same interface is poorly 
defined and appears to be unstable on the simulation 
timescale. One expects the diblock to become homoge- 
neous in the long time limit. 

This preliminary study on diblock copolymers allowed 
us to validate the radical-like copolymerization technics. 
The advantage of this technique resides in the control 
of the geometry of simulated copolymers, as well as the 
possibility to generate in a flexible way configurations 
with various topologies and chain architectures. 

V. CONCLUSION 

1. the radical-like polymer chains generation method 
is inspired by radical polymerisation in which the 
reactive center of a polymer chain consists of a rad- 
ical. The free radical reaction mechanism can be 
divided in to three stages: (i) initiation (creation 
of free radicals); (ii) propagation (construction of 
the repeating chain): (iii) termination (radical is 
no longer active). 

2. Performing an relatively important number of MD 
relaxation steps between each growth step (typi- 
cally 300) leads to well equilibrated chains (in terms 
of gyration radii, Mean Square Internal Distance, 
and Primitive Path Analysis), even for relatively 
short relaxation stage (10 6 MD Steps). 



3. The main advantage of the radical-like generation 
algorithm is that equilibration occurs simultane- 
ously as chain growth, within a coarse grained 
molecular dynamics scheme. 

4. Polymer melts generated with the radical-like algo- 
rithm are as well equilibrated as melts generated 
by more classical methods (like fast push-off). 

5. The radical- like generation method is particularly 
adapted to generate polydisperse polymer melts 
(branched polymers, star polymers, copolymers,...). 

6. nano-structured lamellar block copolymers have 
been successfully generated with the radical-like 
method. 

7. physical and mechanical properties of di-block and 
tri-block copolymers generated using this algo- 
rithm, will be the subject of a futur paper. 
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